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Abstract 

Gamma-ray bursts are the most luminous explosions in the Uni- 
verse, whose origin and mechanism is the focus of intense interest. 
They appear connected to supernova remnants from massive stars or the 
merger of their remnants, and their brightness makes them temporar- 
ily detectable out to the larges distances yet explored in the Universe. 
After pioneering breakthroughs from space and ground experiments, 
their study is entering a new phase with observations from the recently 
launched Fermi satellite, as well as the prospect of detections or limits 
from large neutrino and gravitational wave detectors. The interplay be- 
tween such observations and theoretical models of gamma-ray bursts is 
reviewed, as well as their connections to supernovae and cosmology. 

1 Introduction 

Roughly once a day, somewhere within our Hubble horizon a Gamma-ray burst 
(GRB) occurs, which for the next few seconds or tens of seconds completely 
overwhelms the gamma-ray flux from the rest of the Universe, including the 
Sun. In fact, the GRB prompt electromagnetic energy output during tens of 
seconds is comparable to that of the Sun over ~ few x 10 10 years, or to that of 
our entire Milky Way over a few years; and their X-ray and optical afterglow 
over the first day after the outburst can outshine the brightest quasars, as 
well as supernovae, making them potentially important probes of the distant 
Universe. Since the discovery of their X-ray afterglows by the Beppo-SAX 
satellite in 1997 and the subsequent detection of their optical counterparts, 
we have measured these objects out to the farthest cosmological distances. 
Thanks to triggers and measurements from the Swift pQ and Fermi [2] [3] 
satellites, we have now detailed multi-wavelength data for many hundreds of 
bursts, and redshifts for over 200 of them, and this data set will continue to 
grow with the continuation of Swift and Fermi , and the possible upcoming 
S VOM mission pg. 

GRBs are thought to arise either when a massive star 25M ) undergoes 
core collapse, or possibly when a double neutron star or a neutron star and 
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a black hole binary merges [5]. The first scenario applies to so-called long 
GRBs (LGRBs), whose 7-ray light-curve is lasts for t 7 *5 2 s, while the second 
scenario is the likeliest one so far for short GRBs, whose 7-ray light-curves last 
£y ^ 2 s [6] (for the latter, this short duration refers to photons at e ~ 100 keV; 
some "short" bursts, at softer energies, have tails lasting as much as 100 s 
[HE])- In either scenario, it appears inevitable that a compact core object of a 
few to several solar masses forms, whose radius is of order of the Schwarzschild 
radius for this mass, r g ~ 10 6 (M/3M Q ) cm, over a timescale comparable to 
a few dynamic (free-fall) times, which is likely to be a black hole. Accretion 
of residual infalling gas leads, if the core is fast rotating (guaranteed for a 
binary), to an accretion disk whose inner radius is r ~ 3r g ~ 10 7 cm, and 
the typical variability timescale of accretion is to ~ (2GM/rg)~ 1 / 2 ~ 10~ 3 s. 
The bulk of the gravitational energy, of order a solar rest mass or 10 54 erg is, 
as in SNe, rapidly radiated as thermal neutrinos [E v ^h ~ 10 MeV), and some 
amount is radiated as gravitational waves. A smaller fraction, of order Ej ~ 
10 51 — 10 52 erg is converted into a fireball of equivalent blackbody temperature 
T ~ few MeV. This energy eventually emerges as the burst, in the form 
of a jet. However, for purposes of the dynamics one can generally use (see 
below) the isotropic equivalent energy E = Ej^iir/Qj) ~ 10 53 E 53 erg, which 
with a nominal total burst duration ti ~ 10 s implies a nominal isotropic 
equivalent luminosity L 1 ~ E Q /t b = 10 52 L 52 erg s" 1 (if most of that energy 
is emitted as 7-rays). The number density of photons at r is roughly given 
by L 7 = AnT^cri-yE where e ~ kT ~ MeV, and the "compactness parameter 
(roughly the optical depth of a photon with energy <> m e c 2 against 77 — >■ e + e~ 
pair production) is 

£' ~ r 77 ~ n^a T r ~ ~ 10 15 , (1) 

where o~t is the Thomson cross section and a is the fraction of the luminosity 
above m e c 2 . This creates a fireball of gamma-rays, electron-positron pairs and 
hot baryons, where most of the entropy and pressure is in the photons and 
leptons. The optical depth is huge, and the radiation pressure far exceeds 
gravity, so the fireball expands and becomes relativistic. A simple lower limit 
on the expansion bulk Lorentz factor follows from the observations of photons 
up to ^ GeV energies, in some bursts. Such photons, trying to escape the 
source, would collide against softer photons and pair produce, 7 + 7 — > e + + e~, 
degrading the spectrum to ^ 0.5 MeV. However the pair production threshold 
is angle dependent, and pair production is avoided if 

£i£ 2 (l -cos#) < 2m 2 c 4 . (2) 

In a relativistically moving jet, causality implies that only photons within 
angles 9 & 1/T can interact, so with E\ ~ 30 GeV, e 2 ~ MeV, and cos^ ~ 
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1 — 8 2 /2, we see that this implies 



T £ y/( £l /m e c?)(e 2 /m e c?)/2 £ 10 2 . (3) 

A more general constraint on T is obtained by considering the typical photon 
spectral distribution in GRBs, which is a broken power law "Band" function 
n(e) oc ph cm -3 MeV -1 , where (3 ~ 1 or (3 ~ 2 for e below or above a 
break frequency Eb T ~ MeV [9j [10]. The 77 optical depth at each energy 
> m e c 2 depends on the optical depth to target photons at < m e c 2 satisfying 
the threshold condition (T5]), in the jet comoving frame. This optical depth can 
be shown to be oc T -6 , so for increasingly high T the source becomes optically 
thin to increasingly higher energy photons. The result is that typical GRBs, 
even if the highest energy photons observed are only 100 MeV, require bulk 
Lorentz factors in excess of T ~ 250 [TT] . 

If the entire burst energy is released in impulsively, injecting an energy E in 
a timescale to inside a radius ro, with the numbers comparable to those above, 
the initial entropy per baryon is 77 ~ Eq/Mqc 2 , where Mq is the baryon load of 
the fireball. If the pressure is mainly due to radiation and pairs, the inertia is 
due to baryons ( "baryonic dynamics" regime) the bulk Lorentz factor initially 
accelerates as T(r) ~ (r/ro), e.g. [12]. After the baryons have become non- 
relativistic in their own frame, the expansion changes to a coasting behavior at 
a saturation radius r sat ~ r rj, and the fireball continues to expand freely with 
r ~ 77 ~ constant. The observationally estimated values are r ~ 77 ~ 10 2 — 10 3 , 
so the baryon load is typically 10 -5 — 1O _6 M . The behavior is similar if the 
energy and mass input is spread out over accretion times (i.e. outflow feeding 
or ejection times) of ~ 10 — 100 s, as inferred for "long" GRB. On the other 
hand, if the fireball pressure, or rather stress tensor, is dominated by magnetic 
fields, the dynamic behavior is different; depending on the symmetries of the 
fields, the acceleration behavior can range from Y oc r 1 / 3 , e.g. [131 Q31 [TS] to 
r oc r v where 1/3 < rj ^ 2/3, at least when the outflow is one-dimensional 
[TBI [TT] . This regime is referred to as magnetically dominated, or Poynting 
dominated dynamics. 

In practice, the outflow is inferred to be jet-like, rather than isotropic, 
with an average solid angle (Qj) /4tt ~ 1/500 or (0,) ~ 1/30 [HE]. In the 
case of core collapse ( "long" ) GRBs, this can be due to the outer parts of the 
star providing a massive barrier, which is best pierced along the centrifugally 
lightened rotation axis, along which the fireball escapes. The stellar envelope 
provides a sideways pressure which channels the jet. However, as long as the 
jet opening angle 9j exceeds 1/r ~ 10 -2 (i.e. ~ 0.5°), which is generally the 
case, the expansion occurs as if it were isotropic: causality prevents the gas 
to have any knowledge of what happens outside an angle 1/r. For compact 
binary mergers, the data on jet opening angles is much sparser, but the average 
value may not be too different [19]. 
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2 Observations 



The Swift mission, launched in November 2004, finds bursts and observes the 
prompt phase with the Burst Alert Telescope (BAT). The afterglow is then 
observed with the X-Ray Telescope (XRT) and the UV Optical Telescope 
(UVOT). Measurements of the redshift and studies of host galaxies are typi- 
cally done with large ground-based telescopes which receive immediate alerts 
from the spacecraft when GRBs are detected. Swift has, by far, the largest 
number of well-localized bursts, afterglow observations and redshift determi- 
nations. As of 1 April 2012, BAT has detected 669 GRBs (annual average 
rate of ~ 90 per year). Approximately 80% of the BAT-detected GRBs have 
rapid repointings (the remaining 20% have spacecraft constraints that prevent 
rapid slewing). Of those, virtually all long bursts observed promptly have 
detected X-ray afterglow. Short bursts are more likely to have negligible X- 
ray afterglow, fading rapidly below the XRT sensitivity limit. The fraction of 
rapid-pointing GRBs that have UVOT detection is ~ 35%. Combined with 
ground-based optical observations, about ~ 60% of Swift GRB have optical 
afterglow detection. There are so far (mid 2012) about 200 Swift GRBs with 
redshifts, compared with 41 in the pre-Swift era. The redshift distribution of 
Swift GRBs is shown in Fig. |2j 
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Figure 1: Redshift 
distribution and 
cosmic look-back 
time of GRBs. The 
Swift GRBs are in 
blue, the pre-Swift 
GRBs in yellow 
and the co-moving 
volume of the uni- 
verse is the red 
curve. The GRBs 
roughly follow the 
co-moving volume 



The Fermi mission, launched in June 2008, has two instruments, the Gamma- 
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ray Burst Monitor (GBM) and the Large Area Telescope (LAT). The GBM 
has scintillation detectors and covers the energy range from 8 keV to 40 MeV. 
It measures spectra of GRBs and determines their position to ~ 5° accuracy. 
The LAT is a pair conversion telescope covering the energy range from 20 MeV 
to 300 GeV. It measures spectra of sources and positions them to an accuracy 
of < 1°. The GBM detects GRBs at a rate of ~ 250 per year, of which on 
average 20% are short bursts. The LAT detects bursts at a rate of ~ 8 per 
year. 

3 The photon spectrum: standard picture 

The expansion converts internal energy into bulk kinetic energy, so that the 
gas cools in its own rest frame and soon becomes an inefficient radiator. In the 
absence of dissipation, at the photospheric radius r p h where the flow becomes 
optically thin to scattering (which for baryon dominated dynamics usually 
occurs above the saturation radius r sat ) the escaping radiation would carry 
only a small fraction of the burst kinetic energy, and might be expected to 
have a quasi-black-body spectrum [20| [21] . This motivated the fireball shock 
model, where the bulk kinetic energy is reconverted by shocks into random 
particle energy, and thence into non-thermal radiation, at radii beyond the 
scattering photosphere, where the flow is optically thin; this could be at either 
external shocks [2"2"| |2"3"]. where the jet interacts with external (e.g. interstellar) 
matter, or at internal shocks [23] occurring within the jet, at radii intermediate 
between the photosphere and the external shock radius. These radii can be 
expressed as 

r ph ~ (La T /47rm p cV ) ~ 4 x 10 12 L 7i52 ^ cm, 
r dis ~ T 2 ct v ~ 3 x 10 13 77£ 6 *„_ 2 cm, 

r dec ~ {ZEu/Aim^mpcWY 1 * ~ 2 x W 17 E^n^rff 3 cm. (4) 

Here the photospheric radius assumes baryonic dynamics [21] (for magnetic 
dynamics see [15]). The dissipation (internal shock) radius follows from the 
relativistic relation between observer time t, the radius r and Lorentz factor, 
r ~ T 2 ct. Considering two shells of matter ejected at time intervals comparable 
to the variability timescale of ejection t v , with Lorentz factors differing by 
Ar ~ T, and the deceleration radius r dec where external shocks start follows 
from the energy conservation assumption, E ~ (4-7T / ' 3)r^ ec n ex tm p c 2 T 2 , when 
the swept up matter has been shock-heated to an energy comparable to the 
explosion energy (with r ~ rj, e.g. [TJ] for details). 

In the (collisionless) shocks, the particles are reheated to thermal energies 
comparable to the pre-shock relative kinetic energies per particle. The in- 
ternal shocks are semi-relativistic (since Ar ~ T, the relative Lorentz factor 
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T re i = (l/2)[(ri/r 2 ) + (r 2 /ri)] is semi-relativistic), and this results in a shock 
luminosity L s h = e s hL where L is initial kinetic luminosity L = E f]/tbMc 2 
and the shock dissipation efficiency e s h ~ 0.1. Particles repeatedly bounce 
across the shock, scattered by magnetic irregularities whose energy density can 
be assumed to build up to some degree of equipartition with the proton ther- 
mal energy, so the comoving magnetic field is B' 2 /87r ~ eB4T 2 el m p c 2 (primed 
quantities are in comoving frame), where eg < 1. The repeated crossing Fermi 
accelerates the particles to a relativistic power law [22| [23] ; the minimum elec- 
tron (comoving) random Lorentz factor is 7 ejm ~ e e (m p /m e )T re i ^> 1 (for 
internal shocks T re i ~ 1, while for external T re i ~ rj) and one expects a power 
law iV(7 e ) oc 7~ p above that, with p ~ 2 — 2.3. One also expects something 
similar for the protons in the flow. 

i) Internal shock prompt radiation: The electrons in the internal shock 
will emit synchrotron and inverse Compton (IC) radiation, leading to non- 
thermal broken power law photon spectrum, roughly similar to the observed 
"Band" spectra [231 123] • For reasonable values of L, r), t v , e#, e e the synchrotron 
peak energy (observer frame) corresponding to the minimum 7 e m are compa- 
rable to the observed Band spectral break energies, 

L l/2 

e sy , m ~ e br ~ 1 e B /2 e 3 J 2 2 7 ' 52 MeV. (5) 

V2.5 t v,-2 

Eq. (JSJ) assumes the randomized kinetic luminosity of internal shocks L s h 
to be related to their 7-ray luminosity through L 7 = e e L s , which is true in 
the fast cooling regime where cooling time is shorter than the dynamic time 
t'sync ^ t'dyn ES], which in internal shocks is true [2B]. In this fast cooling 
regime, for a typical Fermi electron index p ~ 2 the photon spectral index 
above Ebr is expected to be n(e) oc e~^ 2 ^ 1 oc e~ 2 ph cm 3 s 1 , as typical 
for the high energy branch of the canonical Band spectrum. The synchrotron 
model predicts below e sy ^ m a low energy branch n(e) oc e~ 2 / 3 , which through 
superposition of maxima for various parameters could fit the observed Band 
average low energy branch n(e) oc s~ x [23J (but flatter spectra are a problem 
see below). 

ii) External shock radiation and afterglow: At the relativistic ejecta 
has used up about half its initial energy in sweeping up an amount M sw ~ 
M e j ecta /r] of external material, driving a forward shock into the external gas 
and a reverse shock into the ejecta. This occurs at an observer time 

W,o ~ r dec /{2cr, 2 ) ~ 10 (^ 53 /n ) 1/3 % 5 8/3 s. (6) 

The forward shock is initially highly relativistic, T s h ~ 77, so T re ij s ~ 7] and the 
synchrotron spectrum is in the hard X-rays or gamma-rays. The reverse shock 
builds up slowly and for usual conditions becomes semi-relativistic, T re [ trs ~ 1 
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at the deceleration time, when it has crossed the ejecta. For this reason its 
7 6j m is smaller than that of the forward shock electrons, and the reverse shock 
spectrum peaks in the optical or UV [271 [22]. Beyond rd ec the expansion 
continues but it is increasingly slowed down due to increasing amount of swept 
up matter. In the adiabatic approximation the bulk Lorentz factor changes 
from being ~ 77 ~ constant to a power law decline behavior given by E oc r 3 r 2 , 
or 

r~r/(r/r dec )- 3 / 2 . (7) 

In both the forward and the reverse shock one expects again Fermi acceleration 
of electrons to a power law distribution leading to synchrotron and inverse 
Compton radiation, but the synchrotron break energy becomes softer in time 
as the Doppler boost decreases in accordance with eq.flTJ). This leads to an 
afterglow [28] progressing from X-rays through optical to radio lasting from 
minutes to days to months, with fluxes decaying as power laws in time. This 
prediction was indeed confirmed by observations with the Beppo-SAX satellite 
of the X-ray [29] and with ground telescopes of the optical [30] afterglow of 
GRB 970228, soon followed by the first confirmation of a cosmological redshift 
[5T] and a radio detection [32] for GRB 970508. The amount of data on, and 
understanding of, afterglows has since increased enormously, see e.g. [7J[S]. 

The external reverse shock gas, most luminous at t ~ t^ ec , is in pressure 
equilibrium with forward shock gas, and having a higher particle density and 
smaller energy per electron than the forward shock, its synchrotron spectrum 
peaks in the O /U V. This was predicted to lead to an observable prompt optical 
emission j27J [2S], later detected with robotic ground telescopes such as ROTSE 
[35] triggered by spacecraft, the number of such detections being now several 
dozen [7J. 

The external forward shock is expected to also give rise to an IC com- 
ponent, in particular a synchrotron-self- Compton (SSC) from upscattering its 
own synchrotron photons [271131], which would appear in the GeV range. Such 
GeV emission was detected already by EGRET [35], and more recently by the 
Fermi LAT, e.g. [35]. This is discussed in §0 

4 Prompt MeV emission: issues and develop- 
ments 

Issues arise with the radiation efficiency of internal shocks, which is small in 
the bolometric sense (5-10%), unless the different shells have widely differing 
Lorentz factors J37J [3H1 EH]- The MeV efficiency is also substantially affected 
by IC losses [101 E], in the BATSE range being typically ~ 1 — 5%, both when 
the MeV break is due to synchrotron [37J 03] and when it is due to inverse 
Compton [H]. 
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The synchrotron interpretation of the GRB radiation is the most attractive; 
however, a number of effects can modify the simple synchrotron spectrum. 
One is that the cooling could be rapid, i.e. when the comoving synchrotron 
cooling time t' = 9mgC 5 / 4e 4 B' 2/ ~f e ) ~ 7x 10 8 /B' 2/ ~f e s is less than the comoving 
dynamic time t' dyn ~ r/2cT, the electrons cool down to 7 C = 6iim e c/ ' aTB' 2 t' dyn 
and the spectrum above v c ~ r(3/87r)(eS'/m e c)7^ is F u oc z/ -1 / 2 [251 EH]- 

The radiative efficiency issue has motivated investigating various alterna- 
tives, e.g. relativistic turbulence in the emission region [4*61 Wf\. This assumes 
that relativistic eddies with Lorentz factors 7 r ~ 10 exist in the comoving frame 
of the bulk Y > 300 flow, and survive to undergo at least 7 r changes over a 
dynamic time, leading both to high variability and better efficiency Various 
constraints may however pose difficulties [48], while numerical simulations [49] 
indicate that relativistic turbulence would lead to shocks and thermalization, 
reducing it to non-relativistic. 

The synchrotron spectral interpretation faces a problem from the observed 
distribution of Band low energy spectral indices /3i (where N e oc e@ l below the 
spectral peak), which has a mean value ~ —1, but for a fraction of bursts 
this slope reaches values > —2/3 which are incompatible with a low energy 
synchrotron asymptote fii = —2/3 [50]. Possible explanations include syn- 
chrotron self-absorption in the X-ray [5T] or in the optical range up-scattered 
to X-rays [44], low-pitch angle scattering or jitter radiation [521 E3], or time- 
dependent acceleration [51] , where low-pitch angle diffusion might also explain 
high energy indices steeper than predicted by isotropic scattering. 

Pair formation can become important [241 HOl [41] in internal shocks or 
dissipation regions occurring at small radii, since a high comoving luminosity 
implies a large comoving compactness parameter £' ^> 1 . Pair-breakdown may 
cause a continuous rather then an abrupt heating and lead to a self-regulating 
moderate optical thickness pair plasma at sub-relativistic temperature, sug- 
gesting a comptonized spectrum [35] . Copious pair formation in internal shocks 
may in fact extend the photosphere beyond the baryonic photosphere value (J4J . 
Generic photosphere plus internal shock models [551 ESI EZ] which includes the 
emission of a thermal photosphere as well as a non-thermal component from 
internal shocks outside of it, subject to pair breakdown, which can produce 
both steep low energy spectra, preferred breaks and a power law at high ener- 
gies. A moderate to high scattering depth can lead to a Compton equilibrium 
which gives spectral peaks in the right energy range [SB]. Pair enrichment of 
the outflow (due to back-scatter 77 interactions) can in general affect both the 
radiative efficiency and the spectrum [591 Ell EU E2l 163] . 
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4.1 Photospheric models 

In the synchrotron interpretation the observed peak frequency is dependent on 
the bulk Lorentz factor, which may be random, and since the observed peaks 
appear to concentrate near 0.2-1 MeV [7J, the question can be posed whether 
this is indeed due to synchrotron, or to some other effect. An alternative is to 
attribute a preferred peak to a black-body at the comoving pair recombination 
temperature in the fireball photosphere [61]. In this case a steep low energy 
spectral slope is due to the Rayleigh- Jeans part of the photosphere, and the 
high energy power law spectra and GeV emission require a separate explana- 
tion [55]. A related explanation has been invoked [65] > considering scattering of 
photospheric photons off MHD turbulence in the coasting portion of the out- 
flow, which up-scatters the adiabatically cooled photons up to the observed 
break energy and forms a power law above. 

For a photosphere occurring at r < r sat , which in a baryon-dominated 
model requires high values of rj, the radiative luminosity in the observer frame 
is undiminished, since E' rad oc r _1 but r oc r so E rad ~ constant, or L ph oc 
r 2 r 2 T' 4 oc constant, since T' oc r _1 . However for the more moderate values 
of rj the photosphere occurs at r > r sat , and whereas the kinetic energy of 
the baryons is constant E^ in ~ E ~ constant the radiation energy drops as 
Erad oc (r/r sat )~ 2/3 , or L ph ~ L (r ph /r sat )~ 2/3 [231 [55]. This weakening of the 
photospheric luminosity leads again to a lowered efficiency, as well as a lower 
peak energy than observed. However, if the photosphere is dissipative (due to 
shocks or other dissipation occurring at or below the photosphere) then a high 
efficiency is regained, and the thermal peak photon energies are in the range of 
observed Band peaks [66]. An important aspect is that Compton equilibrium 
of internal shock electrons or pairs with photospheric photons lead to a high 
radiative efficiency, as well as to spectra with a break at the right preferred 
energy and steep low energy slopes [661 EU [68] . It also leads to possible physical 
explanations for the Amati [69] or Ghirlanda [70] relations between spectral 
peak energy and burst fluence [661 EI] • 

4.2 Magnetic models 

An alternative set of models for the prompt emission assume that this is due to 
magnetic reconnection or dissipation processes, or else to the external shock. 
Magnetic models fall into two categories, one where baryons are absent or dy- 
namically negligible, at least initially [721 [T31 [721 El], and another where the 
baryon load is significant, although dynamically sub-dominant relative to the 
magnetic stresses [651 ESI [63] • These scenarios would in all cases still lead to 
an external shock, whose radius would be again given by r^ec in equ. (jl]), with 
a standard forward blast wave, but possibly a weaker or absent reverse shock 
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[TSJUB], due to the very high Alfven (sound) speed in the ejecta. For the same 
reason, internal shocks may be prevented from forming in magnetized outflows 
However, this depends on the magnetization parameter a\ if not too large, re- 
verse shocks [76], [TTJ ESj or internal shocks might still form [79], although with 
different strengths and radiation characteristics. In fact, "internal" dissipation 
regions may form due to magnetic reconnection, at radii comparable but dif- 
fering from r^is of eq.fllj), where electric fields due to reconnection (instead of a 
Fermi mechanism) leads to particle acceleration, and a high radiative efficiency 
is conceivable. 

A hybrid dissipation model, entitled ICMART [80] involves a hybrid mag- 
netically dominated outflow leading to semi-relativistic turbulent reconnection. 
Here a moderately mag netized a = (B' 2 /Airp'c 2 ) < 100 MHD outflow under- 
goes internal shocks as o — > 1, leading to turbulence and reconnection which 
accelerates electrons at radii r > 10 15 cm. These involve fewer protons than 
usual baryonic models, hence less conspicuous photospheres, and have signifi- 
cant variability, and the efficiency and spectrum are argued to have advantages 
over those in the usual synchrotron internal shock models. 

The baryon-free Poynting jet models resemble pulsar wind models, except 
for being jet-shaped, as in AGN baryon-poor models. The energy requirements 
of GRB (isotropic-equivalent luminosities L 1 > 10 52 erg s" 1 ) require magnetic 
fields at the base in excess of B ~ 10 15 G, which can be produced by shear 
and instabilities in an accreting torus around the black hole (BH). The energy 
source can be either the accretion energy, or via the magnetic coupling between 
the disk and BH, extraction of angular momentum from the latter occurring via 
the Blandford-Znajek mechanism [81]. The stresses in in this type of model are 
initially magnetic, involving also pairs and photons, and just as in purely hydro 
baryon-loaded models they lead to an initial Lorentz factor growth T oc r up to 
a pair annihilation photosphere [82J. This provides a first radiation component, 
typically peaking in the hard X-ray to MeV, with upscattering adding a high 
energy power law. Internal shocks are not expected beyond this photosphere, 
but an external shock provides another IC component, which reaches into the 
GeV-TeV range. 

The baryon-loaded magnetically dominated jets have a different accelera- 
tion dynamics than the baryon-poor magnetic jets or the baryon dominated 
hydrodynamic jets: whereas both the latter accelerate initially as V oc r and 
eventually achieve a coasting Lorentz factor Tf~ L 1 /Mc 2 , the baryon-loaded 
magnetically dominated jets have a variety of possible acceleration behaviors, 
generally less steep than the above. In the simplest treatment of a homoge- 
neous jet with transverse magnetic field which undergoes reconnection, the 
acceleration is T oc r 1//3 [731 US] , while in inhomogeneous jets where the mag- 
netic field and the rest mass varies across the jet the average acceleration 
ranges from T oc r 1 / 3 to various other power laws intermediate between this 
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and r (x r [531 E]- F ew calculations have been made [51] of the expected 
(leptonic) spectral signatures in the simpler magnetized outflow photospheres, 
typically in a one-zone steady state approximation, showing that a Band-type 
spectrum can be reproduced. 



5 GeV-TeV phenomenology and models 

The first Fermi GRB observations, starting in late 2008, soon yielded a number 
of surprises. One of the first bright objects showing radically new features was 
GRB 0809 16C [36], in which the GeV emission started only with a second 
pulse, which was delayed by ~ 4 s relative to the first pulse, which visible only 
in MeV (Fig. EJ. 
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The spectra of GRB 0809 16C consisted of simple Band- type broken power 
laws, the first pulse having a soft high energy index disappearing at GeV, but 
the second and subsequent pulses having harder high energy indices reaching 
well into the GeV range. There was no evidence for a second spectral compo- 
nent (such as expected from inverse Compton or hadronic effects). The peak 
energy of the Band function evolved from soft to hard and back to soft, but 
in this as well as in other Fermi LAT bursts, the GeV emission persisted in 
afterglows typically lasting > 1000 s. On the other hand, in a few bursts, such 
as GRB090902B [85], a second spectral component did indeed appear, at 5a 
significance, and also a lower energy power law extension whose significance is 
lower but suggestive, Another burst with a high energy second component was 
GRB 090926A [86], this one showing a clear cut-off or turnover to the high 
energy power law (Fig. [5]). 
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Figure 3: 
Spectra of 
GRB090926A 
from Fermi at 
four different 
time intervals, 
a= [0.0-3.3s], 
b= [3.3-9. 7s], 
c= [9.7-10.5s], 
d= [10.5-21.6s] 
PI- 



10 10 2 10 3 10" 10 5 10 6 

Energy (keV) 

A significant advance from Fermi LAT was the discovery of the first GeV 
short burst, GRB 090510 [87], whose general behavior (including a GeV delay) 
was qualitatively similar to that of long bursts. Several more short bursts have 
been discovered since with the Fermi LAT. 

The Fermi -LAT extended emission, if one ignores various details, has a 
relatively simple interpretation in terms of conventional forward shock lep- 
tonic synchrotron models (i.e. relying on accelerated electrons or e + e~ pairs) 
[88] [89] . Such models provide a natural delay between an assumed prompt 
MeV emission (assumed implicitly to come from, e.g. internal shocks or other 
"inner" mechanisms) and the GeV emission from the external shock, which 
starts after a few seconds time delay. However, taking into account more care- 
fully the constraints provided by the Swift MeV and X-ray observations, and 
considering carefully the accompanying inverse Compton (IC) scattering and 
Klein-Nishina effects, it is clear that at least during the prompt emission, there 
must be a subtle interplay between the shorter lasting mechanism providing 
the MeV radiation and the mechanism or emission region responsible for the 
bulk of the longer lasting GeV radiation [9U] EU E2]- One general shortcoming 
of these early studies was a postponement of addressing the interaction of the 
GeV emission with a specific, self-consistent model of the prompt emission, 
including the radiative inefficiency in an implicit internal shock assumption. 

A resolution of this problem is possible if the prompt MeV Band spec- 
trum is due to an efficient dissipative photosphere (baryonic, in this case) with 
an internal shock upscattering the MeV photons at a lower efficiency, giving 
the delayed GeV spectrum [93]. Alternatively, for a magnetically dominated 
outflow, where internal shocks may not occur, an efficient dissipative photo- 
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spheric Band spectrum can be up-scattered by the external shock and produce 
the observed delayed GeV spectrum [HI]. Depending on the parameters, the 
combined spectrum can look like a two-component or a single Band spectrum 
(Fig. BJ. On the other hand, a delayed GeV spectrum can also be expected in 
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Figure 4: A magnetically dominated leptonic model where the MeV Band 
spectrum is due to photospheric emission, there are no internal shocks, and the 
external reverse and forward shock upscatter the MeV photospheric spectrum 
into the GeV range. Parameters are typical for Fermi LAT GRBs, but in some 
cases lead (left) to a two-component spectrum, while in others (right panel) it 
can be fitted as a single Band spectrum extending to the GeV range [91]. 



hadronic models, which assume the co-acceleration, along with the electrons, of 
protons which undergo electromagnetic cascades and synchrotron losses along 
with their secondaries [95j EE1 [97J |98] see $61 

The LAT data show that a fraction of GRB are emitting (in their own rest 
frame) photons in the energy range of at least up to 30—90 GeV. A partial list of 
Fermi LAT detections [99] of maximum observed photon energies and redshifts 
(E^ obs , z) is (13.2, 4.35), (7.5, 3.57), (5.3, 0.74), (31.3, 0.90), (33.4, 1.82), (19.6, 2.10), 
(2.8, 0.897), (4.3, 1.37). This list shows that (i) even z > 4 bursts can produce 
E y > 10 GeV photons at the observer, and (ii) some z ~ 1 bursts can pro- 
duce Ey > 30 GeV photons at the observer. This is highly encouraging for 
the planned large Cherenkov Telescope Array (CTA), as described in recent 
reviews [1001 [M]. The CTA detection rate is estimated HOD] to be 0.7 - 1.6 
per year, based on the rate of Swift triggers (while GBM triggers on Fermi are 
more frequent, their positional accuracy is poorer). This rate is affected by 
uncertainties in the fraction of bursts which emit in the GeV range, relative 
to those emitting below 100 MeV [1021 [103]. E.g., as of February 2011, in 
2.5 years, Fermi LAT detected 4 bursts at energies > 10 GeV (or 20 at > 0.1 
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GeV) out of some 700 bursts detected by Fermi GBM at E < 100 MeV. This 
very small fraction of the total (< 1%) of course is in part due to the size 
constraints under which space detectors must operate. 

In the standard internal shock model of prompt emission, the intra-source 
77 absorption typically prevents photons in excess of a few GeV to emerge 
[4"0"| |4"T] , unless the bulk Lorentz factor is above ~ 700 |104] . For photospheric 
models of the prompt emission, e.g. |105j . photons in excess of 10 GeV can 
escape the source from radii r 77 ~ 10 15 cm, and such radii are also inferred 
phenomenologically from one-zone analyses of the Fermi data on GRB. How- 
ever, most of the GeV emission occurs during the afterglow, which is good for 
ground-based TeV Cherenkov telescopes, whose reaction time can be slower. 
Indeed, the GeV emission can last up to ~ 1000 s, far more than the ~ 2 — 50 s 
of the MeV emission. In the standard external shock scenario, the compact- 
ness parameter is smaller than in the internal shock, and inverse Compton 
scattering is expected to lead to multi-GeV and TeV photons [751 EI], the de- 
tails depending on the electron distribution slope and the radiative regime (e.g. 
slow or fast cooling). This scenario is thought to be responsible for the after- 
glows of GRB [28] , and is also thought to be responsible for the extended GeV 
emission observed by LAT so far [881 Ell I106[ 1107] . etc. Of course, propagation 
in the intergalactic medium from high redshifts leads to additional 77 — > e ± 
interaction with the extragalactic background light, or EBL |108[ I109[ 1110] . 
the threshold for which depends on the photon energy and the source redshift. 

Thus, if TeV emission is produced, it is mainly expected to be detectable 
from z < 0.5, while the 10-30 GeV emission should be (and is) detectable from 
higher redshifts. Thus, the GeV detectability is dictated by the source physics, 
the source rate and the immediate source environment. The source rate, based 
on MeV observations, is well constrained [7J, while the near-source environ- 
ment effects can be reasonably parametrized (e.g. |lllj ). The source physics, 
however, has large uncertainties. E.g. in an external shock model the simple 
synchrotron self-Compton (SSC) model can be additionally complicated by the 
scattering of photons arising at other locations well inside the external shock, 
e.g. from the photosphere [03], or from an inner region energized by continued 
central engine activity |112j . Similar uncertainties about the soft photon source 
and location would affect hadronic cascade models. The observed GBM high 
energy spectral slopes are in many cases steep enough not to expect much GeV 
emission from their extrapolation [99], while in other cases the LAT spectrum 
shows a cutoff or turnover, e.g. in GRB 090926B [86]. Nonetheless, all things 
considered, the estimate |100] of 0.7 — 1.6 CTA detections per year appears to 
be a conservative lower limit. 
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6 Hadronic models 



If GRB jets are baryonic, or magnetically dominated but with non-negligible 
baryon load, the charged baryons should be co-accelerated with the electrons 
in any shocks or reconnection zones, and hadronic processes would lead to 
both secondary high energy photons and neutrinos. Monte Carlo codes have 
been developed to model hadronic effects in relativistic flows, including p, 7 
cascades, Bethe-Heitler interactions, etc. E.g. |113[ 196] used such a code to 
calculate the photon spectra from secondary leptons resulting from hadronic 
interactions following proton acceleration in the same shocks that accelerate 
primary electrons in GRBs. The code uses an escape probability formulation 
to compute the emerging spectra in a steady state, and provides a detailed 
quantification of the signatures of hadronic interactions, which can be com- 
pared to those arising from purely leptonic acceleration. Spectral fits of the 
Fermi LAT observations of the short GRB 090510 were modeled by [96] as 
electron synchrotron for the MeV component and photohadronic cascade ra- 
diation for the GeV distinct power law component. 

Since acceleration as well as cascade development can take some time, in 
principle even one-zone models might result in GeV-MeV photon delays. E.g. 
[95] assumes for GRB 090510 the prompt MeV to be electron synchrotron and 
the GeV to be proton synchrotron, whose cooling time cranks down the pho- 
ton energy into the GeV range on a few second delay timescale, the electron 
plus proton synchrotron merging into a single Band function with the approx- 
imate spectral slope of the GeV photons. A more recent one-zone hadronic 
calculation [98] shows that even when proton synchrotron is not important, 
hadronic cascade development leads to a second GeV component, with time 
delays comparable to the observed ones (Fig. El right panel). Similar delays 
can, however, be also obtained in purely leptonic two-zone photosphere plus 
external shock models [114] (Fig. [51 left panel). 

Hadronic interactions can also have implications for a low energy photon 
power law below the Band function, perhaps resulting in a GRB optical prompt 
flash, as discussed by [115] . For the usual Band MeV spectrum produced by 
conventional leptonic mechanisms, the acceleration of hadrons leads to secon- 
daries whose radiation produces both a high energy "extra" GeV component 
and a prompt bright optical emission from secondary synchrotron. This might 
explain, e.g. the observed "naked eye" 5th magnitude flash of GRB 080319B, 
e.g. HE]. 

Hadronic binary collisions in baryon-loaded jets can also be important, 
both for efficient energy dissipation and for shaping the photon spectrum. 
This is because the baryons will consist of both protons (p) and neutrons (n), 
especially if heavy elements are photo-dissociated. The protons are coupled to 
the radiation during the acceleration phase but the neutrons are carried along 
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Figure 5: Temporal evolution of the observable spectral photon flux for typical 
Fermi LAT parameters, from Monte Carlo simulations. Left: a purely leptonic 
two-zone model with a photospheric (MeV) Band component and upscattering 
into the GeV range by a shock further out |114j . Right: a one- zone hadronic 
model, where electron synchrotron produces the Band MeV spectrum and 
hadronic cascade secondaries produce the GeV spectrum, as well as a low 
energy component |98j . 

only thanks to nuclear (p, n) elastic collisions, whose characteristic timescale 
at some point becomes longer than the expansion time. At this point the p 
and n radial relative drift velocity v approaches c, leading to the collisions 
becoming inelastic, p + n — > 7r + ,7r°, in turn leading to positrons, gamma- 
rays and neutrinos |117j . Such inelastic (p, n) collisions can also arise in jets 
where the bulk Lorentz factor is transversely inhomogeneous |118j . e.g. going 
from large to small as the angle increases, as expected intuitively from a jet 
experiencing friction against the surrounding stellar envelope. In such cases, 
the neutrons from the slower, outer jet regions can diffuse into the faster inner 
regions, leading to inelastic (p, n) and (n, n) collisions resulting again in pions. 
An interesting consequence of either radial or tangential (n,p) drifts is that 
the decoupling generally occurs below the scattering photosphere, and the 
resulting positrons and gamma-rays deposit a significant fraction of the relative 
kinetic energy into the flow, reheating it |105j . Internal dissipation below 
the photosphere has been advocated, e.g. |119] to explain the MeV peaks as 
quasi-thermal photospheric peaks |120[ 1121] . while having a large radiative 
efficiency Such internal dissipation is naturally provided by (p, n) decoupling, 
and numerical simulations |1U5] indicate that a Band spectrum and a high 
efficiency is indeed obtained, which remains the case even when the flow is 
magnetized up to Eb = 2 |122j . while keeping the dynamics dominated by the 
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baryons. These numerical results were obtained for nominal cases based on a 
specific radial (n,p) velocity difference, although the phenomenon is generic. 

The photon spectral signatures of a magnetically dominated, baryon loaded 
leptonic plus hadronic GRB model involving nuclear collisions has been cal- 
culated by |123j . This uses a realistic transverse structure of a fast core- 
slow sheath. The analytical results indicate that the transverse neutron col- 
lisions become most effective, resulting in GeV photons at radii from which 
the observer-frame time delay relative to the photospheric MeV photons is ap- 
propriate to explain the observed Fermi time lags. The purely leptonic (SSC, 
EIC) time delays and spectral components of such a baryon-loaded magnetic 
model, in the absence of drifts and transverse gradients, have been calculated 
by |124j . leading to delays in the range observed by Fermi. 

A hadronic model which attempts to self consistently produce the GeV 
radiation, the MeV Band spectrum and the low energy (optical) power law is 
discussed in [125j . The protons accelerated in the shocks or magnetic recon- 
nection regions result in hadronic cascades which produce, as in |115j . the GeV 
and optical power laws, while the cooled leptonic secondaries are re-accelerated 
via Fermi 2nd order mechanism in the MHD turbulent waves produced by the 
same shocks or reconnection regions, leading self-consistently to an MeV Band 
spectrum. 

7 Gravitational waves from GRB 

GRB may be also sources of gravitational waves (GWs). The most likely 
such sources are short GRBs |126j . if these indeed arise from merging compact 
objects [7]. The rates in advanced LIGO and VIRGO may be at least several 
per year [127] . Long GRBs are more speculatively as sources, since in the 
favored core collapse scenario the collapse may be more chaotic [128j . They 
may, nonetheless, be weakly detectable as GW sources, especially if the core 
collapse breaks up into substantial blobs |129j , or if they go through a magnetar 
phase leading to a bar [130J. More recent, detailed numerical calculations of 
collapsar (long) GRBs lead to GW prospects which range from pessimistic 
[T5T] to modest [T32] . 

8 High energy neutrinos from GRBs 

High energy (10 9 eV ^ E u ^ 10 18 eV neutrinos may be expected from baryon- 
loaded GRBs if sufficient protons are co-accelerated in the shocks |133j . The 
most widely considered paradigm involves proton acceleration and p'j inter- 
actions in internal shocks, resulting in prompt ~ 100 TeV HENUs [134[ 1135] . 
Other interaction regions considered are external shocks, with py interactions 
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on reverse shock UV photons leading to EeV HENUs [136j ; and pre-emerging 
or choked jets in collapsars resulting in HENU precursors [137] . Also, for bary- 
onic dominated GRBs, a neutrino component may arise from photospheric pj 
and pp interactions [138[ 1139] . An EeV neutrino flux is also expected from 
external shocks in very massive Pop. Ill magnetically dominated GRBs |140j . 
Current IceCube observations [1414 1142] 1143] are putting significant constraints 
on the internal shock neutrino emission model, with data from the full array 
still to be analyzed. One caveat is that, since the above analysis, several 
groups [144] 1145] 1146] have recalculated the GRB internal shock neutrino pro- 
duction in greater detail, including multi-pion and Kaon production in the 
P7 interactions, and allowing for various astrophysical uncertainties including 
different values of the Lorentz factor and the accelerated proton to electron 
ratio L p /L e = l/f e . The conclusion from these revised calculations is that the 
current IceCube (IC40+IC59) measurements need to be extended for another 
4 to 9 years for obtaining a strong constraint. 

9 Progenitors and Supernova Connection 

Including the collimation correction, the GRB electromagnetic emission is en- 
ergetically quite compatible with an origin in, say, either compact mergers of 
neutron star-neutron star (NS-NS) or black hole-neutron star (BH-NS) binaries 
[2"Ul 1147} 1148} 1149] . or with a core collapse (hypernova or collapsar) model of a 
massive stellar progenitor [150] 1151] 1152] , which would be related to but much 
rarer than core-collapse supernovae 017]. While in both scenarios the outcome 
could be, at least temporarily, a massive fast-rotating ultra-high magnetic field 
neutron star (a magnetar) [1531 [T4] , the mass of the resulting central object ex- 
ceeds substantially the Chandrasekhar mass and is is expected to lead, sooner 
or later, to the formation of a central black hole. The latter will be fed through 
a (seconds to minutes) accretion episode from the surrounding disrupted core 
stellar matter, which provides the energy source for the ejection of relativistic 
matter responsible for the radiation. This inference appears to be confirmed 
by numerical simulations, for NS-NS or NS-BH mergers [154} 11551 1156] as well 
as for collapsar models [152} 1157} 1158] . 

The above numerical simulations also indicate that (1) the compact object 
merger accretion disks are less massive and the accretion episode (when the 
disk is not highly magnetized) lasts less than a few seconds, compatible with 
the observations of the canonical short GRBs (SGRBs); and (2) in the collapsar 
models the accretion lasts for tens of seconds or more, compatible with the 
durations of canonical long GRBs (LGRBs). 

The observations of LGRBs indicate that they are generally located in 
active star-forming environments, usually in blue, small or not too massive, 
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gas-rich galaxies |159[ EJ [7], which is where one expects massive stars to be 
present. Progenitor stars more massive than ~ 25 — 28M , following core 
collapse, are expected to result in a black hole (BH) central remnant, either 
directly or through an intermediate neutron star (NS) phase |1601 1161] . Such 
BH core collapse events, if the core is sufficiently fast rotating, can lead to a 
fall-back fed accretion disk powering a relativistic jet, which is able to escape 
a star of ~ 10 11 cm \157\ 1162] . This radius corresponds to those of Wolf- 
Rayet (WR) stars, which are thought to arise from more massive M > 25M 
progenitors, whittled down by wind mass loss prior to core collapse. The high 
rotation rate, which favors the wind mass loss and also the formation of a 
longer lasting accretion disk, is thought to be enhanced when the star arises in 
a metal-poor environment [5], which in fact seems to characterize many LGRB 
host galaxies |1631 [7] . 

The massive core collapse model of LGRBs is confirmed by the fact that 
LGRBs are, in some cases, demonstrably associated with type Ib/c supernovae, 
whose explosion is, to within errors, contemporaneous with the GRB |1641 
1165} 1166} 1167] . The SNe Ib/Ic are generally thought to have WR progenitors, 
whether they are associated with GRB or not; only a small fraction of order few 
% of SNe Ib/c appear to be associated with LGRBs |168[ I169[ 1167] . However, 
the SNe Ib/c associated with GRBs, as well as a good fraction of those not 
associated, are classified as hypernovae (HNe), meaning that they have unusu- 
ally broad spectral lines indicating a semi-relativistic envelope ejecta (v/c ~ 1), 
and inferred isotropic energies E HN iso ~ 10 52 ' 5 erg, as opposed to the average 
SNe with v/c ~ 0.1 and E SN ,iso ~ 10 51 erg [HQl EH Q221 IT731ITT3] . 

For short gamma-ray bursts (SGRB) the most widely favored candidates 
are mergers of neutron star binaries (DNS) or neutron star-black hole (BH) 
binaries, which lose orbital angular momentum by gravitational wave radiation 
and undergo a merger. This second progenitor scenario has only now begun 
to be tested thanks to the Swift detection of short burst afterglows |175[ I176[ 
I177[ H78[ 1179] . The SGRBs are found both in evolved (elliptical) galaxies |175] 
and in galaxies with star formation [7J 1180] . in proportions compatible with 
that expected for an old population such as neutron stars. While neutron stars 
are expected, and found, in young star-forming galaxies, massive young stars 
are not expected, and not found in old population ellipticals. And, indeed, 
no supernova has ever been found exploding at the same time and location as 
a SGRB [7J. Of course, neutron stars are the product of supernovae, which 
can have occurred much earlier than the burst, from progenitors whose initial 
mass was is 8M & ^ M* ^ 25 M . A kick is imparted to the NS at the 
supernova event, so the NS can wander off significant distances from its birth 
site (many Kpc, even outside the host galaxy). If the NS was born in a binary 
and/or later became a binary, the time between the initial explosion and the 
eventual merger can range up to 10 8 — 10 9 years, and is very unlikely to be 
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less than 10 6 years |181[ 1182] . Only a handful of SGRBs have yielded reliable 
lightcurve breaks suitable for determining a jet opening angle. The latest 
measurements and comparison to previous data [19] indicate an an average 
9j ~ 5° (comparable to the LGRB average value, although there is one outlier 
at 25°). This is interesting, since for DNS mergers there is no stellar envelope 
(as for LGRBs) to provide collimation, at most there would be a wind; however, 
such narrow jets would be compatible with twisted or hoop-stress collimated 
MHD jets from DNS mergers [T56lfT83] . 

10 High redshift GRBs 

Long GRBs are astonishingly bright, both in gamma-rays and at longer wave- 
lenghts. In the optical, typical brightnesses are ~ 18th magnitude a few hours 
after the trigger (and some have been detected in the 5th-10th magnitude range 
seconds after trigger), while a Milky- Way-type galaxy has ~ 32nd magnitude 
at a redshift z — 8. In fact, GRBs vie with galaxies for the record on the 
highest confirmed redshift measurements, e.g. GRB080913 at z = 6.7 |184j, 
GRB090423 at z = 8.2 [H5l[T86] (through spectroscopy), while GRB 090429B 
has a photometric redshift of z ~ 9.4 |187j . It is possible that even much more 
distant objects than these have already been detected in the gamma-ray and 
X-ray detectors of Swift and Fermi, although for such z > 9 objects a specific 
(optical/IR or other) redshift signature is extremely difficult and noisy, so red- 
shift diagnostics are increasingly harder to obtain in this range. The above 
discoveries do, however, indicate that the prospect of eventually reaching into 
the realm of Pop. Ill objects is becoming increasingly realistic. 

Population III GRBs at 10 ^ z ^ 20 may result, as they do at lower 
redshifts, from massive, M > 25 — 3OM metal-poor stars whose core col- 
lapses to a black hole |188j . However, the mass of Pop. Ill stellar progenitors 
could be as high as ~ 1000M Q , leading to 100 - 500M© black holes [185] . 
although the Pop. Ill masses are a subject of debate (and could be much 
lower |190[ I191j . For extremely massive black holes, the jets are likely to be 
Poynting-dominated, e.g. powered by the Blandford-Znajek mechanism. The 
expansion dynamics and the radiation arising from such very massive Poynt- 
ing jet GRBs was discussed by [192] . At typical redshifts z ~ 20 this implies 
a "prompt" emission extending to < 1 day which should be detectable by 
Swift or Fermi , being most prominent initially around 50 keV due to the jet 
pair photosphere, followed after a similar time interval by an external shock 
synchrotron component at a few keV and an inverse Compton component at 
> 70 GeV [193] . Both the 'prompt' emission and the longer-lasting afterglows 
[T93] of such Pop. Ill GRBs should be detectable with the BAT or XRT on 
Swift or the GBM on Fermi. On Swift, image triggers may be the best way to 
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detect them, and some constraints on their rate are provided by radio surveys. 
They are expected to have GeV extensions as well, but redshift determinations 
need to rely on L-band or K-band spectroscopy. 

11 The Cosmology Connection 

Since GRBs are seen out to the largest redshifts yet measured, and for periods 
of hours to days they can outshine any other objects at those distances, their 
potential usefulness as tools for cosmology has been intriguing for some time. 

The simplest way, to use them as distance markers, is unfortunately not 
straightforward. This is because they are not good "standard candles" , which 
could be used e.g. in a Hubble-type diagram to compare flux against redshift 
to compare against cosmological models to deduce a closure parameter, or an 
acceleration rate. Even the collimation-corrected average gamma-ray energy 
i&yj ~ 10 51 ergs has too much variance to be of direct use as a standard 
candle. The error is still almost twice as large as that obtained from SNIa, 
at least at redshifts z & 1.5, which is the most important region for dark 
energy studies. The hope is, still, that one or more of the various empirical 
correlations between observed spectral quantities could lead to a calibration, 
as the Phillips relation does for supernovae, which could turn GRBs into an 
effective distance ruler. 

One such empirical correlation is between the photon spectral peak energy 
E p k and the apparent isotropic energy E iso [69j H94] , namely E p k oc Ef so with 
a ~ 1/2 (Amati relation). Other correlations are between E p k and the peak 
luminosity L p k |195] (Yonetoku relation); or for bursts where the jet opening 
angle is known, between E p k and the collimation corrected gamma-ray energy 
of the jet E^jet |196[ 1197] (Ghirlanda relation); or between E pk and £ 7 j e t 
|198j ; or between E p k, E iso and the light curve break time U r [199] (Liang- 
Zhang relation); or between the X-ray luminosity at break time Lx,br and the 
break time tb r |200j . Of course, these correlations are of interest in themselves 
as possible constraints on the radiation mechanism or emission region, and 
various interpretations have been made, e.g. |2U1} 1119} [7T| I2U2] . However, 
for cosmology only the tightness of the empirical correlation is what matters. 
Observational biases could, of course, pose problems |203[ 1204] , and circularity 
issues may be a concern; the latter, however can in principle be minimized by 
restricting oneself only to directly observed quantities |205] . 

While such correlations make it possible to extend Hubble diagrams up 
to redshifts z & 8 — 9 with GRBs, their usefulness for deriving cosmological 
parameters is still limited. The main reason is that the matter (dark plus 
baryonic) and the vacuum energy densities evolve as p m oc (1 + z) 3 , pv oc 
(1 + z) 3 ( 1+w \ respectively, where w is the equation of state multiplier, p = wp, 
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for vacuum w = —1 being the currently favored value. At present the vacuum 
energy dominates the dynamics, flv,o — (Pv,o/Pcr,o — 0.7 while fl m ,o — 0.3. 
Going back in time or up in redshift, the vacuum energy grows slowly (or not 
at all, if w = —1), while the matter energy grows fast, so that for redshifts 
z > zy = [(£V,o/^m,o)~ 1 ^ 3 "' — 1] ~ 0.5 the vacuum energy becomes negligible 
compared to that of matter. Thus, the information about the vacuum-related 
acceleration and a vacuum equation of state is gleaned mainly from data at 
z ^ 0.5, where the SNIa with redshifts greatly outnumber the GRBs with 
redshift. Some recent papers using GRBs in Hubble diagrams are, e.g. |206] 
and |207j . the latter including 109 GRBs with known redshifts calibrated with 
567 SNIa. While strongly suggestive, and increasingly interesting, the sample 
is still small compared to SNe la, especially at low z, and the dispersion remains 
larger then for SNe, so the usefulness of GRBs as statistical distance indicators 
remains to be seen. 

Hubble time (Gyr) 

Figure 6: Redshift evo- 
lution of the metalicity 
relative to solar values, 
for GRBs shown with 
blue dots and QSOs show 
with open circles. The 
GRB metalicity is on av- 
erage ~5 times larger 
than in QSO. These are 
based on damped Ly- 
man alpha (DLA) spec- 
tral features. The up- 
per horizontal x-axis in- 
dicates the age of the 
Universe (Hubble time) 
[208] . 

GRBs, however, are likely to be unique as beacons for probing the high 
redshift Universe. They are detectable with current gamma-ray, X-ray and 
infra-red detectors out to distances corresponding to the earliest stars for- 
mation epoch |209[ 1210] . and they may provide possible redshift signatures 
|211[ 1212] extending into the 10 ^ z & 20 range. Their strong, featureless 
power law continuum spectrum shining through the intergalactic medium and 
intervening young galaxies or protogalaxies provide a sensitive probe of the 
ionization state, velocity distribution and chemical composition at those red- 
shifts [213[ 1214] 1215] . An example is shown in Fig. 6, indicating the change of 
the metalicity (given by the oxygen to hydrogen abundance ratio) as a func- 
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tion of redshift. These abundances are determined from spectral absorption 
lines in the continuum radiation of GRBs and quasars. The GRB lines come 
mainly from the host galaxy gas in the star forming region where the explosion 
occurred, while the quasar lines arise in random intervening galaxies along the 

ition out to higher red- 
3r metalicities. This is 
forming regions, which 
explosions. 

Figure 7: Cosmic star for- 
mation history |216j . Shown 
are the data compiled in 
[217] (light circles) and from 
Ly-a emitters (LAE) [2T8] . 
Also recent LBG data for 
two UV LF integrations: 
down to 0.2L* (down tri- 
angles, |219j ). and complete 
up to z = 3 (up triangles). 
Swift GRB-inferred rates 
are diamonds, the shaded 
band shows the range re- 
sulting from varying evolu- 
tionary parameters. Also 
shown is the critical p* from 
[220] for C/f esc = 40, 30, 20 
(dashed lines, top to bot- 
tom), see also |221] . 

GRBs also provide an excellent tool for investigating the cosmic star forma- 
tion rate (SFR) of the high redshift universe, and thereby also the rate of large 
scale structure (LSS) formation, out to (so far) redshifts in the 8-10 range. 
This is exemplified in Fig. 7, which shows the star formation rate determined 
through various techniques. The long GRBs are the endpoints of the lives of 
massive stars, and their rate is therefore approximately proportional to the 
star formation rate in general. However, at high redshifts the rates are very 
uncertain, and may be subject to various observational biases. There may also 
be evolutionary biases, such as a dependence of the long GRB formation on 
the metalicity of the host galaxy, which needs to be taken into account, e.g. 
[2131 1215] . 

The most distant GRBs may also provide the only possible probes of the 
era of the first generation of (Population III) stars formed in the Universe 
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[189[ I192[ I193[ 1222] . These relics of the infant universe could be the most 
sensitive probes of the redshift for the start of large scale structure formation, 
with significant implications for the properties of the dark matter. 

12 The future 

Both Swift and Fermi are likely to be functional and return GRB data for 
many years to come. They have orbital lifetimes extending beyond 2025 and 
no critical expendables. The SVOM mission [I] is an approved Chinese- French 
mission to observe GRBs and their afterglow. It has a wide-field instrument 
to image the bursts and one to study the spectrum. The spacecraft rapidly 
slews like Swift to point X-ray and optical telescopes for afterglow observation. 
There are other concepts in consideration, such as Lobster [223], which per- 
forms the wide-field observations in the X-ray band suitable for high-redshift 
GRBs. In addition, combined with such electromagnetic detection programs, 
increasingly sensitive multi-messenger detection attempts will continue to be 
pursued using high energy neutrinos (§8]) and gravitational waves (^7J). These 
expanded observational efforts will require more detailed theoretical interpreta- 
tion and models, extending well beyond what has been achieved so far. Based 
on past experience, GRBs may be counted on to provide further exciting sur- 
prises during the next decade. 

Acknowledgment: We are grateful to NASA for support, and to our col- 
leagues for many useful discussions. 
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